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ÍNDICE GENERAL 


Resumen 


En el presente trabajo se estudia la dispersión de ondas electromagnéticas por dos super- 
ficies esféricas de grafeno concéntricas utilizando un desarrollo basado en las soluciones de 
Mie [Mie08]. Se realizaron cálculos para las secciones eficaces de dispersión y extinción para 
frecuencias de incidencia en el rango de los 'Terahertz donde el grafeno soporta plasmones 
superficiales que se manifiestan como intensificaciones en los espectros de estas secciones 


eficaces. Los datos se verificaron con resultados disponibles en la literatura y se analizaron 





nuevas configuraciones que incluyeron la combinación con materiales dieléctricos y metálicos. 
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RESUMEN 


Capítulo 1 


Introducción 


1.1. Motivación 


El grafeno es un alótropo del carbono en forma de una red cristalina hexagonal de un 
átomo de espesor. Se lo puede considerar como un material bidimensional, especialmente 
para frecuencias en el orden de los Terahertz, donde la longitud de onda en el vacío es 
del orden de 10" nm, mientras que las dimensiones atómicas son del orden del Á, esto es 
1071 nm. El premio Nobel de física del año 2010 fue otorgado a Andre Geim y Konstantin 
Novoselov por el descrubrimiento del grafeno (año 2004) y su posterior caracterización. 
Hasta el momento de su descubrimiento la existencia de láminas de carbono de pocos átomos 
de espesor y estables a temperatura ambiente no se creía posible. Los antecedentes cono- 
cidos eran los fullerenos, estos consisten en superficies de carbono cerradas o semicerradas 
de un átomo de espesor, de forma esférica o cilíndrica. El primero de ellos, descubierto en 
1985, es llamado Buckminsterfullerene y consiste en una superficie icosaédrica de carbono 
de un átomo de espesor compuesta por unos 60 átomos [HWKS85| y por cuyo descubri- 
miento recibieron el Premio Nobel de química en 1996 Harold W. Kroto, Robert F. Curl y 
Richard E. Smalley ¡Rea96]|. Luego se han conseguido estructuras esferoides de 70, 76, 78 y 


84 átomos de carbono. En 1993 Sumio lijima descubrió los nanotubos de carbono [ST93]. 





Estos nanotubos también pertenecen a la familia de los fullerenos, en este caso se trata de 
estructuras de un átomo de espesor con forma de cilindros (Figura 1.1). Debido a que estas 


estructuras contienen pocos átomos (del orden de 100), el estudio de su respuesta óptica no 


2 CAPÍTULO 1. INTRODUCCIÓN 


admite un tratamiento macroscópico. En cambio, el grafeno en las frecuencias consideradas 


en este trabajo, se puede tratar de manera macroscópica a través de parámetros constituti- 





vos a determinar en forma experimental o teórica. 





Figura 1.1: Fullerenos: a) Buckminsterfullerene b) nanotubo de carbono 





Algunas de las propiedades que hacen del grafeno un material interesante desde el punto 
de vista tecnológico y teórico son su transparencia (tiene una transmisividad cercana a la 
unidad) y su conductividad. La conductividad del grafeno se puede controlar a través de 
un campo eléctrico transversal. El mecanismo, llamado Efecto de Campo Eléctrico, también 
está presente en los metales, aunque debido al apantallamiento (en distancias del orden de 
1 nm), el efecto sobre la densidad de portadores queda minimizado por las cargas en volumen 


y el cambio en la densidad de portadores resulta menor que el 1%, aún para láminas de 





muy pequeño espesor [KGN04]. El Efecto de Campo Eléctrico se puede conseguir en grafeno 
dispuesto sobre un sustrato dieléctrico mediante una diferencia de potencial entre ellos 
llamado Potencial de Compuerta (Gate Voltage). Se ha observado experimentalmente que 
de esta manera se consigue una densidad de portadores que varía en forma lineal con la 
diferencia de potencial aplicada [KGNO05]. Esto hace que la conductividad del grafeno pueda 


ser ajustada mediante la aplicación de una diferencia de potencial (Figura 1.2). 





Estas propiedades han generado en los últimos años una carrera en el estudio teórico y de 


las posiblidades tecnológicas del grafeno [Dep16|]. 
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MK 


Figura 1.2: Conductividad del grafeno en función de la diferencia de potencial aplicada (Gate 
Voltage) para grafeno sobre $102, tomado de [KG.NO5]| 


1.2. Propiedades constitutivas del grafeno 


Los átomos de carbono en el grafeno se disponen formando hexágonos regulares. Este ti- 
po de red, llamado panal de abeja, no constituye en sí mismo una red de Bravais, pero puede 
describirse mediante una red de Bravais triangular con una base [AM”76|. Para describir la 


ven 


red triangular, se pueden escoger por ejemplo los vectores aj = a(5,v3) y 43 = a(—-5, 4 








en este caso el parámetro de red a es 1,44 A [EVO6]|, y un vector que representa el segundo 


átomo de la base, b = a(0, 75). (Figura 1.3) 


La red recíproca es una red hexagonal, rotada en 30% con respecto a la red directa. Los 





vectores que generan la red recíproca se pueden escribir bi = (1,7) y bo = Z (1, — 7), 





resulta conveniente escribir además los vectores que unen un átomo de la red con sus tres 
primeros vecinos, estos vienen dados por ci = a(0, 73), Ca = a(>3,- 375) ME a(=>5,- 373): 
Utilizando un esquema de cálculo de Tight-Binding [AM76|, se llega, para las energías de 


las bandas electrónicas a [FV06]: 


e1.2(k) = +11/1 + 4cos(3k,a/2)cos(k,yav3/2) + 4cos?(k, av 3/2) (1.1) 
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Figura 1.3: Red de grafeno formada por átomos de carbono con los vectores de la red a; y 
as y el vector de la base b 


donde y = 3eV. 


La primera zona de Brillouin representa también un hexágono en la red recíproca, en la figu- 





ra 1.4 se pueden ver las bandas de energía sobre esta zona, donde se destacan los puntos K en 
los cuales las bandas de valencia y de conducción forman conos que se tocan en sus vértices. 


Debido a esto, el grafeno puede ser considerado un semiconductor con banda prohibida nula. 


En los puntos K de la primera zona de Brillouin, donde la energía es cero, se puede hacer 
un desarrollo para vectores de onda cercanos, considerando en la red recíproca un vector 
k = K + q, donde |q]| es pequeño. El hamiltoniano en términos de q resulta en una relación 


de dispersión lineal para pequeños desplazamientos de los puntos K. 


020) =+5 alg (1.2) 
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Figura 1.4: Bandas de energía del grafeno, el hexágono representa la primera zona de Bri- 





llouin donde se pueden ver los puntos K, en estos puntos se tocan las bandas de valencia y 
de conducción 





Esta relación de dispersión lineal en el vector de onda q posee la forma funcional de la 
dispersión de fermiones sin masa con una velocidad v = 3ya/2 = 10% cm/s en lugar de la 


velocidad de la luz. 


Para frecuencias altas, mucho mayores que kv, la conductividad del grafeno viene dada 





por la siguiente expresión conocida como Fórmula de Kubo [Fal08]: 


_ we le] dfole) e fol=e)— fole) 
Sá ea 1 (Fis)? de de / E ER Pe [es 








donde fol(e) es la función de distribución de Fermi: 


1 


exp [(e — Ap)/KgT] + 1 e 


fole) = 


en esta última expresión, y se refiere al potencial químico y T' a la temperatura absoluta. 


La primera integral en (1.3) representa la conductividad debida a las transiciones dentro 
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de una misma banda de energía o intrabanda, integrando se obtiene: 


dj In[2cosh(£5)] 
E A 1. 
pa huso + ¿hr pla) 


gMtra en 


en la que y representa la frecuencia de colisiones de los electrones con impurezas en la 
superficie y Kg a la constante de Boltzmann. 

La segunda integral en (1.3) representa la conductividad debida a las transiciones entre 
bandas o interbanda [Fal08| 


l 1f1 1 hw — 2hp 0 (hw + 2hp)* 
inter == se — t íIOAA — 1 —_ —_ === 1. 
d (o) 4 (3 j ió ( 2KgT ) ÓN a — 2hp)? + (2KgT)? 20) 


La conductividad del grafeno es la suma de ambas contribuciones: 
o(w) =ca [a ""*(y) + a (y) do) 


donde c es la velocidad de la luz en el vacío y aparece la constante de estructura fina 


EE E 
a = 137: 


A bajas temperaturas, el término intrabanda (1.5) domina la conductividad como puede 


verse en la figura 1.5. 


hue 
2KBgT >? 


hue 
2KgT 








A muy bajas temperaturas, el factor In|[2cosh( )| se puede aproximar por en 
ese caso, el término intrabanda que domina la conductividad del grafeno en estas condiciones 


resulta: 
UE 2 Hule 


1 
E sl 


que tiene la forma de la conductividad clásica derivada del modelo de Drude. 


1.3. Plasmónica y grafeno 


Los plasmones superficiales son ondas en la densidad de carga, al igual que en otros 
fenómenos ondulatorios, estas se pueden presentar como ondas propagantes, en cuyo caso 
se habla de plasmones propagantes o como ondas estacionarias, en cuyo caso hablamos de 
plasmones localizados. 


Los plasmones propagantes aparecen como modos a lo largo de guías de ondas y en 





nanotubos o estructuras cilíndricas y se pueden propagar a lo largo del eje o en espirales 
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Conductividad del grafeno 


Im(0) 











Ll L Ll L L L L 
0 1000 2000 3000 4000 5000 0 1000 2000 3000 4000 5000 


Temperatura (K) Temperatura (K) 


Figura 1.5: Contribución de los términos inter e intrabanda a la conductividad del grafeno 
en función de la temperatura para una frecuencia fija de 1'T'Hz, con un potencial químico 
hu. =1,0eV y frecuencia de colisiones hy. = 0,0001 eV. Como se puede ver, en estos rangos 
de temperatura, tanto la parte real como la imaginaria de la conductividad del grafeno son 
dominadas por la conductividad intrabanda. 





alrededor de la superficie cilíndrica. En los últimos años se han estudiado este tipo de 
plasmones es nanotubos metálicos [SRO8| y en cilindros recubiertos con grafeno [YGJ14a, 
YGJ14b]. Estos plasmones propagantes, a diferencia de los plasmones localizados, tienen un 
espectro de frecuencias continuo. 

Los plasmones localizados, al igual que las ondas estacionarias en muchos sistemas ondu- 
latorios, tienen un espectro de frecuencias discreto [WLK15, Mai07]. En estructuras cilíndri- 
cas y nanotubos, se pueden excitar ambos tipos de plasmones, propagantes y localizados, 
dependiendo de la dirección de incidencia y la frecuencia [MARD15, MCD16). 

Uno de los desafíos actuales de la óptica es controlar la luz en volúmenes de dimensiones 
mucho menores que la longitud de onda al cubo. La plasmónica, que ya ha jugado un papel 
muy importante en la construcción de metamateriales y en el desarrollo de sus aplicaciones, 
se presenta como unos de los caminos más prometedores para tal fin. Los metales nobles 
como la plata o el oro han sido hasta el momento los materiales plasmónicos más usados 
en numerosas aplicaciones. Sin embargo, debido a la imposibilidad de controlar la densidad 


electrónica de los metales nobles, los plasmones metálicos tienen aplicación en un rango 
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acotado del espectro electromagnético. Por ejemplo, no tienen aplicación en la zona de los 





terahertz, donde se esperan novedosas aplicaciones [Sak05| en áreas como medicina, segu- 
ridad o astronomía submilimétrica. La posibilidad de generalizar al rango de los terahertz 
muchos dispositivos basados en la excitación de plasmones superficiales, actualmente en uso 
en el rango visible o infrarojo, ha despertado el interés por el estudio de los plasmones 
superficiales existentes en una lámina de grafeno [MJB13]. 

La estructura electrónica de la capa monoatómica de grafeno tiene zonas donde la relación 
de dispersión es lineal y el gap entre la banda de valencia y la de conducción es nulo 
[KGNO5, ZQLBOS, Fal08]. La densidad de carga del grafeno (y en consecuencia también sus 
propiedades ópticas [XLN13, dA14]) pueden modularse mediante dopaje químico [HLZ11| 
o mediante la aplicación de un campo estático como se mencionó al comienzo de la sección 
[CFCW11, FWS08]. Por estos motivos, el grafeno parece ser un material especialmente 
adecuado para manipular la luz eléctricamente, tal como sucede en numerosas aplicaciones 
donde la excitación de plasmones superficiales juega un rol fundamental [Y VBP10, WGX13.. 

Estudios recientes [MJS09] han mostrado que si se bloquean las transiciones interban- 
da, los plasmones de una hoja plana de grafeno pueden estar muy confinados y a la vez 
propagarse con muy baja pérdida. El bloqueo de las transiciones interbanda puede lograrse 
empujando el umbral de estas transiciones a frecuencias mucho mayores que la frecuencia 


de trabajo, por ejemplo aumentando el nivel de dopaje. 


1.4. Desarrollo de este trabajo 


En este trabajo estudiamos plasmones localizados en partículas esféricas con un núcleo 
metálico o dieléctrico con distintas configuraciones de recubrimiento que involucran una o 
dos capas de grafeno y una corteza dieléctrica. Partiendo de un desarrollo basado en la 
soluciones de Mie se establecen las condiciones de contorno que deben cumplir los campos 


eléctrico y magnético en las distintas interfases. Estas condiciones vinculan a los campos 








en las diferentes regiones del espacio y permiten calcular los coeficientes de transmisión 





(capítulo 2). Estos coeficientes de transmisión relacionan a la onda incidente con el campo 
dispersado y a partir de ellos se derivan diferentes cantidades de interés como las secciones 
eficaces o eficiencias de dispersión, absorción y extinción. Se desarrolla un programa de 


computadora para calcular numericamente estas cantidades, se verifica la convergencia de 
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los cálculos y el correcto funcionamiento (capítulo 3). Finalmente se hallan cantidades como 
las secciones eficaces variando distintos parámetros y se obtienen distribuciones de campos 
cercanos utilizando configuraciones que no han sido estudiadas con anterioridad (capítulo 


4) y se extraen conclusiones de lo realizado (capítulo 5). 
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Capítulo 2 


Desarrollo Teórico 


Damos un resumen de las ecuaciones de Maxwell, que gobiernan el electromagnetismo 
macroscópico en medios lineales, isótropos y homogéneos. La superficie de grafeno será tra- 
tada como una hoja de corriente a la que caracterizaremos por una conductividad superficial 
o(w), que tiene contribuciones intrabanda e interbanda, según el modelo de Kubo que ya fue 
explicado en el capitulo anterior (1.5), (1.6), (1.7). Luego usaremos las ecuaciones de Max- 


well y las condiciones de contorno para analizar partculas del tipo ncleo-recubrimiento (core- 





shell) en las que aparecen capas esféricas de grafeno entre diferentes materiales, dieléctricos 
y metálicos. Las partículas estudiadas se esquematizan en la figura 2.1 donde €, €2 y €3 
representan las permitividades dieléctricas en las zonas r < Ry, Ri <r < Ra y r > Ro 
respectivamente en las que se divide el espacio y R¡ y Ra representan los radios de las capas 
de grafeno. La onda incidente es una onda plana cuyo vector de onda apunta en la dirección 


z positiva y el campo eléctrico tiene polarización lineal en la dirección 2. 


2.1. Fcuaciones de Maxwell 


Se parte de las ecuaciones de Maxwell en el sistema gaussiano: 


V-D=d4rp (2 


V-B=0 (2.2) 
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Figura 2.1: Configuración estudiada. En la figura se pueden distinguir tres regiones del 
espacio. Cada una de estas regiones es homogénea y posée una permitividad eléctrica propia, 
estas son £1, €2 y €3. Dividiendo el espacio en las tres regiones se encuentran las capas 
esféricas de grafeno de radios Ri y Ra que poseén ambas las mismas propiedades, descriptas 
por una conductividad superficial d. 


Ar 10D 
H = — Si — Ze 
V x le 5 (2.3) 
109B 
He = 2 2.4 
Ds c Ot (2.4) 
con 

B =yH (2.5) 


D =«E (2.6) 
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J=0E (247 


donde se considera un medio isótropo y homogéneo, esto es, tanto e como y son escalares 





que, dentro de un mismo medio, no dependen de la posición. De aquí en más ej se referirá 


a la permitividad del medio j. Lo mismo para las permeabilidades 1; 


Si los campos tienen una dependencia temporal armónica, del tipo exp (—iwt), entonces 


las ecuaciones (2.1), (2.2), (2.3) y (2.4) pueden escribirse: 


V.H=0 (2.8) 
V-E=0 (2.9) 
dr WE 
VxH=" 0E-¡—E (2.10) 
C C 
VxE=:“H (2.11) 
C 





en las que ahora E y H pueden ser vectores complejos y no contienen ya la dependencia 
temporal, se dice que el conjunto de ecuaciones anterior representa las ecuaciones de Maxwell 
en el dominio de frecuencias, cuando no hay cargas libres. En general, e, y y la conductividad 
del medio y son funciones de la frecuencia ww. 


En un medio no conductor, a = 0, combinando las ecuaciones (2.10) y (2.11) se tiene: 








2 
Vx (VxE)= E (2.12) 

y 2 
Vx (Vx H)= — (2.13) 


de las que usando la identidad vectorial V x(Wx A) = V(V-A) —V-(VA) y las ecuaciones 


(2.8) y (2.9) resultan las ecuaciones de onda o ecuaciones de Helmholtz homogéneas: 
VE +k“E=0 (2.14) 


V"H +k%“H=0 (2.15) 
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2 . . . . 
con k? = eu que deben satisfacer los campos en ausencia de cargas y corrientes libres. 


2.2. Soluciones de las ecuaciones de los campos en coor- 


denadas esféricas 


Debido a que estamos interesados en estudiar configuraciones con simetría esférica, re- 
sulta conveniente utilizar coordenadas esféricas con vectores unitarios definidos como en la 


(Figura 2.2). 


> 
“> 


Figura 2.2: Versores del sistema de coordenadas que vamos a utilizar a lo largo de este 
trabajo. 


Se puede construir una solución para las ecuaciones (2.14) y (2.15) que cumpla además 


con (2.10) y (2.11) a partir de [Dor10]|]: 
N =-—rx Vy (2.16) 


donde 4 es una función escalar y el signo menos se ha puesto por conveniencia [BH98]. 
Del cálculo vectorial: 


—=rxVy=V x (ry) —- y V xr (2) 


El último término del lado derecho de la ecuación anterior es nulo, ya que el vector posición 
r es irrotacional. Además, al ser N un rotor, su divergencia es nula. 


Introduciendo N en la ecuación de Helmholtz, y utilizando que V x (V*ry) = V x (rV%p), 
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se tiene: 


VÍN + kN = V x [r(V%p + k*9p)] (2.18) 


Entonces, N es solución de la ecuación de Helmholtz vectorial si 4 es una solución de la 
ecuación de Helmholtz escalar. En coordenadas esféricas: 

1 9 ( ,01 1 8) 0v A 5 

—= =— —=—=TI 0) = —=—h.=5+ki=0 2.19 

Or (> 7) y r2sen(0) 06 (sent 3) iS r2sen(0) Op? A) Aa 
Si buscamos una solución en variables separadas, de la forma y = R(r)0(0)9(6) que separe 


la ecuación en: 


d2$ 





qa FP =0 (2.20) 
de a (seno) 7) she e eS 00 O:=U (2.21) 
E (e) + [kór* — L(0+ 1)R =0 (0:00) 


Las soluciones de estas ecuaciones son bien conocidas y figuran en gran parte de la literatura 
sobre electromagnetismo y física cuantica [BH98]. 
Se obtiene para P((): 

$ = exp (img) (2.23) 


La condición de que Y sea univaluada restringe los posibles valores de m a números reales 
enteros. 


Pidiendo que las soluciones de (2.21) sean finitas en 0 = 0 y 0 = 711, las funciones O(0) son 





combinaciones lineales de las funciones asociadas de Legendre de primera especie Pj”(cos(0)) 


con £ = |[m|,|m|+1,|m|+2,.... Estas funciones tienen las propiedades de ortogonalidad 
. 2 ((+m) 
e Er = ==> —, Oo 2.24 
'B A E AT Em e (2.24) 
1 / 
P;"(cos(0))P;” (cos(0)) (L4+m)! 
—————— HA 0)) = ————, Ómm! az 
/ sen?(0) 10)) m(l — m)! pp 


La solución de la parte angular de la ecuación (2.19), Y/”(0,9) = A0(0)9($), donde A es 
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un factor de normalización, está dada entonces por los armónicos esféricos: 


(2£ + 1) (£ —- m)! 





OS ar Cos(0)) exp (imó) (2.26) 
que cumplen con la ortonormalidad: 
277 TE ; 
/ / Y¿" (0,9) Yer” “(0,p) sen(0) dd de = Omm 010 (2.27) 
o JO 
y la relación de completitud: 
00 L 
YY Y7"(0,9) Y¿"(0",9") = ó[cos(0) — cos(0)J0($ — $”) (2.28) 
£=0 m=-—£ 


En la ecuación radial (2.22), se puede hacer la substitución R(r) = r7!uy(r), de manera 
que queda: 


3 +2 +k- 


dr2 - rdr y2 (2:29) 


p 1d A ar) =0 


que representa una ecuación de Bessel con argumento semientero para uy(kr). Las soluciones 





para la parte radial vienen entonces dadas por las funciones de Bessel esféricas de kr 


jelr) = 4/5 Jw+ 4 (kr) (2.30) 


junto con la funciones de Neumann 


yelkr) = 4/5, Yes 3 (er) (2.31) 


y las funciones de Hankel 


pg (ur) = af az [des 3 (br) + Yo, 3 (or) (2.32) 


son combinaciones lineales de las de Bessel y Neumann y representan ondas esféricas salientes 


(+) y entrantes (—). Si x = kr, las funciones esféricas de Bessel y Neumann se pueden obtener 


sea) = al (5 ay a) (2.33) 


a partir de: 





x dx Xx 
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ve) = a (5 a (22) (2.34) 


Entonces, las funciones N¿,, = —r Xx V[fo(kr)Y/"(0,9) = felkr)[—r x VY¿"(0, 6)], donde 


f¿ representa a las funciones de Bessel, Neumann o Hankel esféricas, son soluciones de la 





ecuación de Helmholtz vectorial (2.14) o (2.15). 





Utilizando un procedimiento similar, se pueden definir las funciones My», = V x Nm que 
también verifican las ecuaciones (2.14) o (2.15). 


Se definen los armónicos esféricos vectoriales: 


X2m(0,6) ES A X VYem(0, $) (2.35) 


de manera que en términos de éstos y factores de normalización y conveniencia, las funciones 


N y M quedan definidas como: 








Nem(r) =— Jeri em (F) (2.36) 
Mim(r) = o V X Nom(r) (2.37) 
N 5 (1) = hg (kr) X em(P) (2.38) 
M7, (1) = o V x Nz,, (1) (2.39) 
N?,, (1) = ye(kr)Xem(P) (2.40) 
M2, (1) = ÉV x Ni, (r) (2.41) 
En componentes: 
o Ema a 22+1(4-m)! [mP;"(cos(0)) - ¡IP (cos(0)) » 
Nem(r) = Qe+1) felz) dr (L+m) ( sano qn de ó) 
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Ximo oa) [00 + 1)] PE (cos(0))+ 
((+m)! dd 


Ha DE + costó] (243 


—1 ca el m 


SE 
s 
SS 
eL 
a 
3 
Ele 
j 








con x= kr y fe(x) representando a las funciones de Bessel o de Hankel, según corresponda. 


Algunas propiedades de los armónicos esféricos vectoriales que utilizaremos más adelante 





son: 
dr 
/ (Fx Xogm):Xg1,,, dO =0 (2.45) 
AT 
AT 
Además, se puede ver que: 
AT 


para todo £, L' y m, m' y por su construcción, las N,,, sólo tienen componentes tangenciales, 


mientras que las M,,, tienen tanto componentes tangenciales como componente radial. 





Dada la completitud de las funciones de Bessel esféricas y los armónicos esféricos, cualquier 


solución para la ecuación de Helmholtz escalar puede ser desarrollada como una suma infinita 





de combinaciones de estas funciones y cualquier solución de las ecuaciones de los campos 





puede ser desarrollada como una suma infinita de las funciones Nm y Mm. En la figura 
2.3 mostramos los vectores N(0,¿) y M(0,¿H) sobre la superficie de una esfera. Nótese que 


sólo los vectores M tienen componentes radiales. 
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Figura 2.3: Algunas de las funciones N¿ m y Mo m sobre la superficie de una esfera, sólo las 
funciones M tienen componentes radiales. En la versión electrónica se recomienda ampliar 
la figura para verla con mayor detalle. 


2.3. Desarrollo de los campos 


1. Los campos en el espacio 3, exterior a las dos esferas, son las sumas de los campos 


incidente y dispersado. 


Ez == Pine En sp (2.48) 


Ha = H:ne + Adisp (2.49) 
2. Para los campos incidentes utilizamos una onda plana. 


E(r) = Eo expli(k . r — wt)|] (2.50) 
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A partir de la Ley de Faraday-Lenz: 


p OH 
VxE=-=-=— ZO 
c Ot ) 
. € € 
H=-¿—VxE=>— kx E (2.52) 
pu pu 
como |k| = ./peZ, resulta la siguiente relación entre el campo eléctrico y el campo 


magnético: 


H= EkxE (2.53) 
pl 


. Los campos incidente y dispersado resultan: 


Enval = 2 ón Mm (Y ) + rd Nem(r ) (2.54) 


Hinc(1) =. e Ds Nom(r) — a Mp (1) (2.55) 
Edisp(r 20 Ms (0) + Agp Nm (E) (2.56) 


ES 
Haisp(Y nes dl N; 6 JS epi Mn (1) (2.57) 


Las formas (2.54) a (2.57) para los campos eléctrico y magnético fueron utilizadas por 
primera vez por Gustav Mie en 1908 [Mie08] y el tratamiento de la dispersión por 


una esfera que se utiliza a lo largo de este trabajo se conoce como dispersión de Mie. 





Los coeficientes ad% y alA en los desarrollos (2.54) y (2.55) se eligen de manera que 





los campos incidentes representen una onda plana. Este desarrollo se muestra en la 


próxima sección. Los campos dispersados contienen funciones de Hankel esféricas de 








primera especie hj (kr) que representan ondas salientes. 


. En el espacio 2, entre las dos capas esféricas, los campos contienen las funciones esféri- 





cas de Bessel de primera y segunda especie je(kr) e ye(kr) en su dependencia radial. 


= 2 Mem (1) + 03% Nom(5) + ag72 Mi, (1) + ag00 N7,, (1) (2.58) 
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== o Nom(r) — 3% Mon (1) + YE NY, (1) — 2H M2, (1) (2.59) 
12 





5. Finalmente, en el espacio 1, interior a la capa esférica de radio R¡, la dependencia radial 








viene dada por funciones de Bessel esféricas de primera especie je(kr) solamente, ya 





que las de segunda especie y¿(kr) divergen en el origen 


= 20m Mem(r) + a34 Ne (1) (2.60) 


la E ya E No (1) — al Mya, (1) (2.61) 


de esta manera, las formas que vamos a utilizar para los campos en cada subespacio de 
nuestro arreglo vienen dadas por las ecuaciones (2.54), (2.55), (2.56) y (2.57) para el campo 
en el espacio exterior a las esferas, que denominamos 3. Para el espacio entre las esferas, 
que llamamos 2, los campos vienen dados por las ecuaciones (2.58) y (2.59) y finalmente, 


para el espacio dentro de la esfera interior, que denominamos 1, por (2.60) y (2.61). En 





la subsección siguiente hallaremos los coeficientes de las ecuaciones (2.54), (2.55) para que 





representen una onda plana incidente. 


2.4. Desarrollo de una onda plana en armónicos esféril- 


cos vectoriales 


Al haber simetría esférica, la dirección y polarización de la onda plana incidente se 





pueden elegir de manera arbitraria sin perder generalidad. Entonces vamos a considerar una 
onda que se propaga en la dirección Z, con polarización en la dirección X. En coordenadas 


cartesianas, esto es: 


Eine = Ep e" s0) z (2.62) 
y el campo magnético se obtiene a partir de (2.53). En coordenadas esféricas: 


A 


£ =sen(0) cos(9) + +cos(0) cos(9) O — seníd) $ (2.63) 
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En términos de los armónicos esféricos vectoriales, los campos incidentes vienen dados por 


las ecuaciones (2.54) y (2.55). Si multiplicamos escalarmente cada una de estas ecuaciones 


* 


?m) USamos la propiedad de ortogonalidad con los M.,,, e integramos en un ángulo 


por N 


sólido de 47, obtenemos de la ecuación (2.54), para el coeficiente o 


/ Eine Ny d9 
Un = (2.64) 
Nom: Nz,,, dO 
AT 
En el denominador se tiene: 
No», Ni, / o Ol (2.65) 
AT AT 


donde se ha usado (2.44). 
Para calcular la integral en el numerador de (2.64), podemos escribir los armónicos esféricos 


vectoriales X ¿,,, en componentes a partir de (2.35). 


e O tdci 

KM mM mcos(8)) gim A — ¡ Jm  dAP¿cos(0)) pmo db (2.66) 
L(2 + 1)sen(0) L(£+1) d6 

con 


(2£ + 1)(£ — m)! 


dr (L4+m) (200 


Lem NN 


la integral en el numerador de (2.64) resulta entonces: 


Eo Lem je(kr) / - cos(0) A — Pi(cos(9)) eme — 
L(£ +1) Ar sen 
dP"(cos(0)) 


— ¿senío) 20 


sd eErcost) sen(0) d0do (2.68) 


Al integrar en O, por la ortogonalidad de las funciones trigonométricas, tenemos solamente 


m <= +1, de manera que resulta 


. TT p+l ] 
aL + cos(9) P*(cos(0)) + O sent0)| etr0do” (2:69) 
AÑ 0 
Esto es: 
En Í likr)r f* d 
EU 79 [ESento) PE (cos(0))] eFrooas ii 
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Para la expresión (2.70) que corresponde a m = +1, se puede usar que |AS65|: 


Pa) = ("A aPyr? E [Pe) (2.71) 


donde P¿(x) es el polinomio de Legendre de grado £, en nuestro caso tenemos 





P, Y 
P(cos(0)) = ereoate)) (2.72) 
de 
Los polinomios de Legendre satisfacen la ecuación: 
d dP; 
»7 (sento) ) = —L(£+1) sen(0) Polcos(0)) (2.73) 
(ecuación (2.21) con m = 0). Entonces, (2.70) para m = +1 resulta: 
Eo Lo jelkr) 7 'N nicostó) 
MAAMMMM L£(£+ 1) sen[0) Prlcos[(0)) e *" “MO E! 
df, UC+1) sentó) Pacos(9) (2.4 
que utilizando la siguiente relación [BH98] : 
/ sen(0) Pe(cos(0)) er reos) go = 25 jo (kr) (2.75) 
0 





se puede escribir, luego de sustituir /7, en la expresión (2.74): 


Eov?2t+1 2 7 : ; 
LLE Ger) je (kr) 
ae ES 1 PE ES TE Eo ¿e y 2 +1 VT (2.76) 


jelkr)je(kr) 
El coeficiente correspondiente a m = —1 se calcula utilizando que: 
a (£- m)! 
O RO 2 LS 
O (2.77) 





obteniendose la misma expresión para ese coeficiente. 
ay = EU VUE 1 yr (2.78) 
Para calcular el coeficiente a?% podemos usar la ecuación (2.55). 


0H. 


Procediendo de manera análoga a lo que se hizo para a,,, : 
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Hine = €3 Eo e kr cos(0) y (2.79) 
p3 
con 
Y = sen(0) senló) * +cos(0) senld) Ó+cos[4) $ (2.80) 


y en este caso se obtiene: 


a = FE 47 yY2+1 yr (2.81) 


De manera que los campos correspondientes a la onda plana incidente resultan: 


Ein. = Eo yr y 224 16 (Mo 1 — My 1) + N¿-1 + N2 1] (2.82) 
£=1 


Hao = Eo Vr [2 y VIA a Noa) (Mea Ma] 2 


Recapitulando, nuestro problema se reduce al de encontrar una relación entre los coe- 


: oa +E _+H 1E ,1H _,2E ,2H 
ficientes de la onda incidente de y ae y los coeficientes aj, , Gp, ) Ugo) Apio Uno Um 


2y E 2yH : 
Ap, Y Ap, que describen los campos en las diferentes regiones a través de las expresiones 


(2.56), (2.57), (2.58), (2.59) (2.60) y (2.61). 





2.5. Condiciones de Contorno en las Interfases 


Para las condiciones de contorno que deben cumplir los campos en las superficies de 
las capas esféricas, voy utilizar las ecuaciones de Maxwell en su forma integral para medios 


isótropos y homogéneos: 


dp <E-5da=4m /// pav (2.84) 


dp ¡H-ñda=0 (2.85) 
5 


d E-de==3, ff non do (2.86) 
€ dt JJs 

A 1 
hna To |) ss b nas+2 /f nomas (2.87) 
E C Ss cdt Ss 
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Siendo gd la conductividad superficial y óg una delta de Kronecker, que vale uno sobre la 


superficie de la interfase y cero fuera de ella. 





mt Sint Mint 


Figura 2.4: Diagrama de la superficie conductora separando dos medios interno y externo 
de diferente permitividad eléctrica y permeabilidad magnética. En el esquema se puede ver 
el circuito utilizado para obtener las relaciones entre los campos en los diferentes espacios. 








Vamos a aplicar las ecuaciones (2.84) y (2.85) al volumen en forma de pastilla de la figura 
2.4. El volumen de esta pastilla se puede hacer tender a cero tomando el límite cuando su 
altura tiende a cero. De esta manera, la integral en volumen del lado derecho de (2.84) se 


transforma en el area de la superficie por la densidad superficial de carge ps. Entonces, se 





obtienen las siguientes dos condiciones para las componentes perpendiculares de los campos 


E y H en la interfase. 


O A DA (2.88) 
y 
Mt O (2.89) 


En la que E, y HA, se refieren a las componentes de los campos perpendiculares a la interfase. 


A continuación aplicamos las ecuaciones (2.86) y (2.87) al circuito de la figura 2.4. 
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Tomando nuevamente el límite, ahora cuando la altura A del circuito tiende a cero, el área 
encerrada en el circuito tiende también a cero, con lo que la integral del lado derecho de 


(2.86) resultará nula. Esto implica la continuidad de las componentes tangenciales del campo 





eléctrico en la superficie de la esfera, y puede escribirse: 


Fx (Box — Eint) =0 (2.90) 





En la que hemos vuelto al sistema de coordenadas esférico donde F es el vector unitario 


radial, perpendicular a la superficie. 


De (2.87), tomando t un versor en la superficie, mediante el mismo procedimiento obte- 


nemos: 


á 4 á 
Ax P) > (Hess — Hips) Al = Ho Ena 4 AL (2.91) 
6 


donde dada la continuidad de las componentes tangenciales del campo eléctrico, se ha es- 
cogido el campo tangencial interior para representar esta componente del campo sobre la 
superficie. 


Luego de reordenar los términos resulta: 


4 a 
Px (Her — Hins) — —o Eúns | +8 =0 (2.92) 
C 
La ecuación anterior implica: 
a AT 
PX (Ho: 2 H:;n+) — — 0 Ent (2.93) 
C 


donde la componente tangencial del campo eléctrico viene dada por Ey = E — (E - F)F. 





Estas son condiciones generales que deben cumplir los campos a ambos lados de una inter- 
fase con conductividad entre dos medios, en las siguientes subsecciones aplicaremos estas 


condiciones a los campos en las regiones con permitividades e,, €2 y €3 que están separadas 





por cáscaras esféricas de grafeno para encontrar una relación entre los coeficientes de la onda 


e 0H OE : edo : 1E 1H 
incidente a,,, y Gp, y los coeficientes de los campos en las distintas regiones, a;,, Y Up, 





2H ,2yE 


+E 
Lm > Un, 


O a Pa en la región 3. Esto nos 


e, 2 H ER 
en la región 1, a72, a y as, en la región 2 y a 





permitirá tener expresiones para los campos en las tres regiones a partir del campo debido 


a la onda plana incidente. 
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2.5.1. Condiciones de contorno en la capa esférica interior r = RR; 


e Condición de contorno (2.90) en la superficie esférica interior, de radio [R;: 


ds Po (EV dollar Km) + a jalar) Xen 
+ ao Ñ x (o Xx volar) Xen) al a P Xx yelkor)X gm — 
2 


- [Esizs Je Vx jeta) Xom) + Ai PX ¿0 =0 (2.94) 
Utilizando que 


(2.95) 


x (Vx felkriXgnm) = — (e sd da) X om 


donde fo(kr) representa las funciones je(kr), ye(kr) o hi (kr) según corresponda y 
multiplicando escalarmente por X;,,,, e integrando en un ángulo sólido de 47 y utili- 
zando las propiedades (2.44) y (2.45), sólo sobreviven los términos correspondientes a 


multipolos magnéticos que contienen los coeficientes a ea , y de ,, asociados a 





un modo transverso magnético TM. Se puede entonces despejar la relación entre los 


coeficientes: 


1 jelkaR1)  [Oje(k2r) Ll o2yn [ yelk2R1) , | Oyelkor) 
ko 5 on (Feo Ri il | Or r=R; me ko “em Ri ve Or MR 


lo 15 [JelKkiRa) , | Oje(kar) 
E = Ze 
ra Ue hi E ea Lts 0 ( 96) 


Si definimos x;¿ = k¡R;, entonces la ecuación anterior puede escribirse 


211 [221 jelxo1) ar +01 [291 ye (291) aj — 221 [011 joel) ai =0 (2.97) 


Volviendo a la ecuación (2.94), ahora siguiendo un procedimiento similar multiplicamos 


por F x X;,,,,, utilizamos que 
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en la que el último término es nulo, dado que las Xy,» sólo tienen componentes tan- 





genciales. Si integramos nuevamente en el ángulo sólido e invocamos las propiedades 
(2.44) y (2.45), sobreviven ahora sólo los términos que contienen a los coeficientes 


yH . / . ., 
do da a. , y rd , asociados a un modo transverso eléctrico TE. Ahora la relación 





entre los coeficientes resulta: 


jelca) ae + yelza1) aye, —jele11) ase =0 (2.99) 


Condición (2.93) en la superficie esférica interior, de radio R;: 


E jelkar)  Ojplkar) 
Xom + — === | Xp 
9) aria kar) Fx Xgm + pa En ( p E > tm T 


k Oy (k 
+ 079, yal kar) Fx X em + E yelhar) + Oyelkor) X tm | — 
ko r Or 


fer E Jelfir) , OjeQar) Ñ 
a 2 atea kr) Fx X¿m + a Al a ay X 2m 


' ; (kr 
E Puri x je kar) X em + as5 del kr)X gm — (ass L(L + DA a) 











Donde se ha utilizado que 





P.E=->) amv LU +1) o (2.101) 
L.m 


en la que fp representa a je, y O hi según corresponda. 
Nuevamente utilizamos el procedimiento de multiplicar escalarmente por X;,,,, e in- 


tegrar en un ángulo sólido de 47. Los términos que contienen F x Xy,, desaparecen 





debido a las propiedades de perpendicularidad ya que los Xy,, no tienen componentes 
radiales. 
Para ver que el término que contiene V x je(kar)X¿,, se anula al multiplicar e integrar, 


debemos utilizar primero que: 


Vx felkr)X em = dEl A o A felkr) (1 V"Yen — VYim) (2.102) 


UF) 
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Al multiplicar por X;,,,, e integrar, el primer término desaparece, al igual que el 
término que continene rV*Y,,,, ya que es radial. Utilizando la forma explícita de los 
armónicos esféricos vectoriales X7,, en términos de los armónicos esféricos Y,,, dada 
por (2.35) para el término que contiene el producto con VY¿,», obviando constantes 


de normalización, se tiene: 
Va (ER Y pra Y ro Vi) (2.103) 


El gradiente de los armónicos esféricos es: 


LO 9 A 1 E 
Vo: = N En exp (img)0 + a E PURO exp imojó) 
(2.104) 
donde N es un factor de normalización y las P;” son las funciones asociadas de Le- 


gendre. 


Para la integral de (2.103) en el ángulo sólido tenemos entonces: 


24 1 2 [¿(m — md] md po (cos(0))+ 
L m 0 


OP," (cos(0)) 


+m PI (cos(0)) 5 


| sen(0)d0do (2.105) 


Al integrar en Q sólo quedan los términos tales que mm = m” de manera que se tiene: 








y trim! [7 ARE CcosOp27 cesó) y y 20800 [o op] 
2 0 





y? 00 ON 
(2.106) 
con x = cos(0). 
Los polinomios asociados de Legendre, P/”(x), valen cero en x =-—1 y x =1 cuando 


m 0, por lo tanto la expresión anterior se anula en ese caso, y como m?” multiplica 
a la suma, también vale cero cuando m/ = 0, por lo tanto, toda la suma resulta nula. 


Entonces, al multiplicar escalarmente (2.100) por X;,, en integrar en el ángulo sólido 
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de 47 queda: 


¿E . 
VE EN [221 jelo1) lar + 011 [221 yla Vaz] 
2 


/ez . Ar . l 
= o [211 jeloi)l az AR a 91 111 Taj jole11) erp =Ú (2.107) 
1 


Partimos nuevamente de (2.100), pero ahora multiplicamos escalarmente por F+xX5,,,., 


e integramos en el ángulo sólido de 47. Utilizando la propiedad (2.45), los términos 


2yH 
Lm 


1H 


Da son los que desaparecen esta vez. El último término de 


que contienen a irá y a 


(2.117) desaparece debido a que estamos multiplicando por una función que sólo tiene 
componentes tangenciales. Entonces queda además de los términos que contienen a 


F Xx X¿m, el término que contiene de la forma (2.102). Aquí tenemos: 


Fx X ¿m = —=—=f0((r) (TV Yen 7 VY tm) 


EE (PX) (2.108) 


Ofe(kr) 1 
| Or 





Del producto con el primer término resulta O fp(kr)/Or luego de integrar. El producto 
con el término que contiene rV?*Y,,,, es nulo, ya que es radial y estamos multiplicando 


por un vector tangencial. Finalmente, del último término, tenemos: 


[440 PY) (FxX, yd (2.109) 


El ass 


Aquí vamos a usar que: 





A AGN 2.110 
¿ a ee 
ds 
VYém  VYi id = sn  ibiic (2.111) 
r 


0 


Entonces, la expresión correspondiente a (2.108) resulta: 


Oje( kar) Y joel kar) 


Or 2 = Lo joa (2.112) 
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Luego, resulta la relación entre los coeficientes: 


e , 
a EN je(2a1) e + 311 yeluo1) có a 
V pa 


E 4. 
Y Eds jelt11) ag — 40 [ejela)l Ag =0 (2.113) 
1 


2.5.2. Condiciones de contorno en la capa esférica exterior r = R> 


e Condición de contorno (2.90) en la superficie esférica exterior, de radio Ra: 


20m PE V x jes) Xsm) +ayló Px jelkgr)Xom+ 
+ e Pe Fox (E Xx hi Car) Xom) + e P ox hi (kgr) X em — 
= [20m e YX joker) Xom) Haro Xx jelkor)Xim+ 
Haro $ x (o x Uber) Xen) +as" xy kar)Xem| =0 (2.114) 


Repitiendo el procedimiento utilizado para la capa interior, al multiplicar por X5,,,,, € 
integrar se obtiene para los coeficientes: 


xa2 [132 jolzz2)l'ap + 222 [232 hi (232) az, — 


— 232 [192 jelzo2)l aja — 032 [122 yo(292)l'az, =0 (2.115) 


si se multiplica (2.114) por + x X;,,, y se integra, se obtiene: 


je(u32) aq + hi (032) ajo — ela) añ — yela22) ago =0 (2.116) 


32 CAPÍTULO 2. DESARROLLO TEÓRICO 


e Condición de contorno (2.93) en la superficie esférica exterior, de radio Ra: 


le Ñ / k3r Ojo (kgr 
ez Y [esa kgr) FOX X£m + ra a E Da +) X om + 


¡ ho (kgr)  OhF(k 
es ho (kgr) Px Xem + E e e + a) Xen - 


E jelkar)  Ojplkar) 
9) Xx pa e E o 
aia kar) Fx X¿m + po En ( a SE a tm+t 


o (Y), Qin) y] 





Ñ / 
+as2 ye(kor) PxXXon +4 


> 7 r Or 
Ar Ll 9p ' SE 1 2yE 
E 70 Un YX Jal K2r Ximo + Aim jel Kar) X tm + yA YX Yako) Xen + 
2 2 


lm 


ok k 
+ yo (kar)Xem — (ads an DA EE DY) 3 =0 


(2.117) 


Si de nuevo multiplicamos escalarmente por X;,, , e integramos en un ángulo sólido 


de 47, se tiene: 


e : 
(E losa less 020 ca lesa 142) 
e , 
ON ES 199 jelto2)) az + 232 [129 yola) aj a 
2 


AT. 
+ E 09 Lo 132 je(222) A + ye (122) a =0 (2.118) 


Si ahora multiplicamos por F x X;,, , en integramos: 


m 


E 
4/ E 292 Je(L32) ALE + 29 hi (232) al = 


E2 : 
+ ma ES je(xo2) rd + 222 Ye(L22) ce => 
V pa 


4. 
= a 02 [22 je(u22)l tm + [122 yel222))' cel =0- (2.119) 


Esto completa los dos conjuntos de cuatro ecuaciones para los coeficientes de los campos 








en términos de los coeficientes del campo incidente 0 y as El primer conjunto relaciona 





los coeficientes que contienen un supraíndice E, que se puede demostrar que corresponden 


a ondas en un modo transverso magnético y el segundo,con supraíndices H a un modo 
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transverso eléctrico de los campos de radiación [BH98, Dorl10]. 


2.6. Coeficientes de Transmisión 


2.6.1. Modo transverso magnético 


Las ecuaciones que relacionan los coeficientes correspondientes a cada modo, vienen 





dadas por (2.97), (2.113), (2.115) y (2.119). Para encontrar el coeficiente de transmisión T/* 


que relaciona los coeficientes de la onda plana incidente y la onda dispersada en el modo 





transverso magnético, debemos resolver el sistema de ecuaciones: 


211 [291 jelxa1)l ar + 211 [221 yelxo1)]/ajyo — 221 [011 jelx11)] aj =0 (2.120) 


e ¿ 
Ue EN je(x21) da + u11 yelxo1) cel = 
2 


le E dr. 
o un je(211) Am — Fi [rrjelor1)) az, =0 (2.121) 
1 


xo [132 jelzz2)]'apy + 222 [132 hi (232) a, — 


— 232 [2022 joltos) din — 232 [222 yola) apo, =0 (2.122) 


E y 
Vi [129 jelzz2) aj + 222 hi (232) aj, | — 


E j 
= Lo ES Je(%22) ie + 299 Ye (u22) | E 
V pa 


Ar . 
PE [o jelx22)1' ajo + [1992 ye( 229) seed =0 (2.123) 


En lo que sigue voy a utilizar la notación die; para referirme a las expresiones |x;¿j0(2;5))”, 
dy;, para [uijyel(xi5)! y dh;, para [u;¡h] (x;5)]/. Las dos primeras, a menos de un signo, 


representan las derivadas de las funciones de Ricatti-Bessel [AS65]. 


yl = 
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2y E 
tm * 


md 








A partir de la ecuación (2.120) se puede despejar a;,, en términos de e y a 


E 
e xiidjó,ad£ + e11dyf, as”. (2.124) 


2210jÍ7 


Reemplazando la expresión anterior en (2.121) se obtiene: 


E ] 
E (mu jo021)az 3 onylon)azy) == 
V pa 


-£ .2E £ _2yE 
El 2110j31 UG, Al 111dY310p>, 
= y del) AAA A A AA 
al 


2310, 





2y 
2110j31 Gio + U11dY51 0 


AT E 
— —o011d)j 2.125 
e 114711 2a1djf, | ( ) 





yE 


2 
Esta última ecuación puede resolverse para a), 


en términos de a definiendo la relación 
entre estos coeficientes como iaa = A o resulta: 


lares )ardj + mjor) adj + Foridji, 11052, 


A A O 
my a1)221dj, — a - Toridjf a11dyS; 


Utilizando que 211 = ,/€1M1k0R1 = ,/€111%01 y una relación análoga para x21, dividiendo el 





numerador y el denominador por 211%01 y usando que 2q1 = ¿R¡. Definiendo g¡ (w) = a : 
resulta: 

H _ —esjelaar)dji, + erjel11) dj + 9 (w)djj1 dijo (2.197) 

mo coyelaa dit, — erjel 11 )dyl, — djt, dyl 

oYelx221)dj1, — Erjelc11)dys, — g1[w)djidysy 
Reemplazando ad n (2.122) y despejando a%,, se tiene: 
dif OE dhf +E 
ad _ 122 113920 Pon, =p 99 320 Pin, (2.128) 


232dj57 + tzadysy AR 


Llm 


Reemplazando en (2.123), definiendo el coeficiente de transmisión T' E de manera que Pol = 





> an se tiene para el coeficiente: 


V 1 U22U32 dj3> + dy32 A) jelt32) — V e Udo jelx22) + yelo22) Ain] digo — Eozito [digo + dia Ain] digo 


a 1 U22032 dio + dy3s Abi ] dh 239 je (22) En ye( 222) Af dhsz + O gÍoo dio Ay dys2 Af] dh) 


(2.129) 


Utilizando como antes que 2;¿= y/E¡M5koHR¿, dividiendo el numerador y el denominador 
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de la ecuación anterior por 2p2222 resulta: 


€3 [dj39 + dys2 Ain | jelt32) — €2 [jolta2) + yolta2) Aza] digo — ga (w) [dizo + dy32 Afin | di3o 
—€3 ldj32 + dyo Ají, ) dh3o + 2 [jolx22) + ye[222) Ap, | dhza + ga[w) [diaz + dy3o Ag, | do 
(2.130) 


H" 
a 


2.6.2. Modo transverso eléctrico 


Por otro lado queda un segundo conjunto de cuatro ecuaciones, (2.99), (2.107), (2.116) y 
(2.118) para el modo transverso eléctrico (los coeficientes que llevan el supraíndice H). Las 


ecuaciones a utilizar en este modo son: 


jelca) ají + yelza1) ao, —jele11) ale =0 (2.131) 


€2 -£ 2H £ _2yH 
1/ 6 EN dj310 gm + T11 lY21 Ue, == 
2 


el : Ar. Ñ 
5 a diia ón ES a 01 211 221 Je(X11) da =0- (2.132) 
1 


jelu32) az + hi (232) al — jelo2) az — yelz22) ajo, =0 (2.133) 


E3 ] 
4/ sn [229 djz20g,, + 222 dh32a7,, | — 
3 
E . H 
24] a ES diégarA + aso dyógapo. | 28 


AT . : 
+ an 09 To 132 jee) > + ye (122) -redl] =0 (2.134) 


de (2.131) resulta: 
sa dela + ye(221)apo, 


E : 2.135 
Llm de (211) ( ) 
Reemplazando la expresión anterior en (2.132), despejando esa ecuación para ad en térmi- 
nos de ds definiendo As, de manera que a = AS ea , resulta ahora: 
1; O a E : d 
A 7 eL 11 dior — y Jel21)dji1 + 91) jelar)je(r1) (2.136) 


Je (11) dy, =F e Ye(221)dj1, — g9i(w)jel221)je(211) 


T 


¡O yl o O RN ME A TE E A A de e A o A A 
“— —Ldhbo [jelm22) + yela2 AP] + Lhi (232) [dio + dyó¿ AF] — gal) jah (232) [jela22) + yela22) AF! 
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Substituyendo en (2.133): 


om Jelts2 aa + ha, (32)0 4, 


a m . ZAS 
y ju(x292) + yelzo2) Aj ) 


Finalmente, substituyendo en (2.134), se tiene: 


sa Yss [jelz22) + yelaz2) AF] — io jelxa2) [djzo + dy3s AF] + gol) ro je lazo [joloo2) + yelz22) A 





En esta sección hemos obtenido los coeficientes de transmisión T/F ec. (2.138) y TF, 


+H 


em <on los co- 





ec. (2.130) que relacionan a los coeficientes de la onda dispersada ed y a 


rrespondientes a la onda incidente de y de Estos coeficientes de transmisión acoplan 





los multipolos, representados por £, de la onda incidente a los multipolos del mismo orden 
en la onda dispersada y no dependen del índice m. Todas las cantidades de interés como 
secciones eficaces, campos cercanos, campos de radiación, etc, se determinan a partir de los 


coeficientes calculados. 


2.7. Flujo de Energía 


Nos interesa conocer los flujos de energía a partir de la onda incidente, esto es, cuánta 
de la energía incidente es irradiada en las ondas salientes y qué porción es disipada en cada 
capa esférica debido al efecto Joule. 

El promedio temporal del flujo de energía a través de una unidad de superficie viene 
dado por: 

E 


(S) = ¿RefE x H') (2.139) 


T 


Luego, la potencia que atraviesa una superficie esférica de radio r se obtiene integrando 
este vector en la dirección perpendicular a esta superficie en cada punto, sobre la superficie 
esférica: 


P= | RAExH')-?r2d0 (2.140) 
ST 4r 





Se puede definir una potencia de extinción, de aquí en más P..,., como la potencia que es 


removida de la radiación incidente, entonces: 


Pext = Pabs + Paisp (2.141) 


(2.138) 
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Utilizando las expresiones de los campos dispersados (2.56) y (2.57) en (2.140), se puede 
obener la potencia dispersada como el flujo del vector de Poyting debido a los campos 
dispersados, saliente a través de una superficie esférica de radio r mayor que los radios de 
las esferas de grafeno. Luego de hacer una permutación en los productos escalar y vectorial, 


resulta: 


Ps — ] Re((Fx Edisp) * Hip) 12d0 (2.142) 


Utilizando propiedades de los armónicos esféricos vectoriales, el producto vectorial resulta: 


¡ pik Oh (k 
is», cat E + a) Xóm +05 hi (kr) Fx Xem| (2.143) 
L m 


Y para el complejo conjugado del campo magnético: 


Hiisp = V ua E a hy (kr) Xn + 


Oh; (de 
e e, X Ktm + hy (er) (1V?Y go, — VY Zn) | (2.144) 





UFO 


donde se ha utilizado (2.102) y que hF*(kr) = h; (kr). 

En virtud de (2.44), (2.45) y (2.98), las integrales que contienen productos entre los armóni- 
cos esféricos vectoriales X ¿y y/o F x X¿2,. son sencillas. Por otro lado, las integrales que 
contienen el producto con LY son cero, ya que los armónicos esféricos vectoriales X ,,,, 
y Fx X£,y no tienen componente radial. 

Las restantes integrales tienen las formas (2.103) y (2.109). 


Resulta entonces: 


Or 


+ PER y Ez + O) hi (er) pe (2,145) 


c Je ho (kr Ohi (kr _ 
Paisp = 3 yy Pe E SIT al E + ma) hy (kr) + 


donde quedan productos de la forma hj(kr)h; (kr) que resulta real, por lo tanto no van 
a aportar a la potencia dispersada y los productos de las funciones de Hankel de primera 


especie con la derivada de las funciones de Hankel de segunda especie y vice versa. Tomando 
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el primero de ellos y escribiendolo en términos de las funciones je(kr) e y¿(kr) resulta: 


Ohj (kr) 


a ia deter) Ger) Oye (kr) a) 


SI O lt) (tr) — a 
(2.146) 





Donde sólo la parte imaginaria aporta a la potencia dispersada. Esta parte puede rescribirse: 


A CUR) 4, (o) — e E) =ikljelkr)ys(kr) —Goer)yelter)] (2.147) 


la cantidad entre corchetes es el Wronskiano de j¿ y ye, usando que Wlje(p),yelp) =p"? y 


haciendo lo mismo con el otro producto, se tiene: 


Paso = oz A] rela (2.148) 


De manera similar se puede calcular la potencia absorbida, como el flujo entrante en la 





misma esfera de radio r, del vector de Poynting debido a los campos totales en el espacio 
exterior a las superficies de grafeno [Jac99]. Esto es: 


Pr E Re (Eme + Edisp) < (Ho + Hiso))) ? 1240 (2.149) 
T 


inc 
4 


(Jue utilizando una permutación cíclica de los productos de vectores, resulta: 


¡AN == Re[[? x (Eme + Exisp)] - (Hipo + Hijisp)) 1240 (2.150) 
AT 


El producto de los campos incidentes resulta imaginario puro, de manera que no aporta a 
la potencia absorbida. Los otros productos que aparecen en la expresión (2.150) son, por un 
lado el producto de los campos dispersados, cuya integral resulta igual a (2.148) pero ahora 
con un signo negativo, ya que corresponde a la potencia entrante en la superficie, y por otro 
lado los productos cruzados entre campos incidentes y dispersados. En los cálculos aparece 
de nuevo el Wronskiano entre las j¿ y las y¿ de Bessel con el que se simplifican algunas 


expresiones para obtener: 


Ex Hx E H 
a ba dE Pelada +aen n ye laz lr + laz el (2.151) 
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Luego, utilizando (2.141), se tiene para la potencia exinguida: 
Po=-£ 817 ros yr y 0H pH 2.152 
ext — J:2 y et 04m Ue Un Ue J ( : ) 


La potencia que lleva la onda plana incidente (2.62) y (2.79) se puede calcular a partir 


del vector de Poynting de esta onda, y resulta: 


S=,/ El (2.153) 
43 


2.7.1. Energía disipada en las capas de grafeno 


En cada una de las capas de grafeno se disipa energía por efecto Joule. La potencia 


disipada en cada capa viene dada por: 
PE 1 J- E RÍdO = / OE? RÍdQ (2.154) 
dr AT 


Donde R; se refiere al radio de la capa 1 = 1,2. En esta cuenta entra la parte del campo 





eléctrico tangencial a la superficie de la esfera E. 
Tenemos: 


[E | = ¡7 X E| (2.155) 


Donde el lado derecho tiene la forma (2.143), de manera que para la potencia disipada en 
la superficie de una capa esférica obtenemos: 
1 E 
Prep =D pega VERO FAROS E labs? + F206RO. E labs, (2.156) 
L m 
Hasta aquí nos hemos concentrado en estudiar de forma teórica rigurosa la dispersión 


dada por dos esferas concéntricas de grafeno sobre las que incide una onda plana. Á partir 


de aquí presentaremos algunos resultados numéricos de interés. 
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Capítulo 3 


Herramienta numérica 


Para la evaluación numérica de los campos, secciones eficaces, etc. se desarrolló un pro- 
grama en lenguaje Python. Con el objetivo de verificar el correcto funcionamiento del pro- 
grama se reprodujeron resultados hallados en la bibliografía sobre el tema. [|XF'S14, TCMI5, 
BYZ15]. Todos los programas utilizados para el desarrollo del código cumplen con los re- 


querimientos de software libre [Hun07]. 


Verificamos el correcto funcionamiento del código y la convergencia de los cálculos com- 
parando los resultados obtenidos con algunos hallados en la literatura, estos resultados 


involucran un núcleo dieléctrico y con un recubrimiento de grafeno. 


3.1. Verificación del funcionamiento del código 


e Esfera dieléctrica sin recubrimiento con distintos valores de índice de refracción [XF'S14]: 
La primera verificación se llevó a cabo considerando una esfera dieléctrica de permi- 


tividad uniforme sin grafeno y variando el índice de refracción de ésta con los valores 





dados en [XFS14|. En la figura (3.1) se pueden ver los resultados obtenidos para las 
secciones eficaces de dispersión divididas por la sección eficaz geométrica, que es igual 
a 217. R” en función de un parámetro de tamaño igual al número de onda por el radio 
de la esfera kyR. En este caso se dejó fijo el radio de la esfera en 10nm y se varió la 


frecuencia de la onda incidente. Los valores obtenidos están en perfecto acuerdo con 
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los dados en la bibliografía. Se puede ver en la figura (3.1) (c), que cuando el índice 
de refracción se vuelve muy grande, la sección eficaz tiende a un valor que es el doble 


de la sección eficaz geométrica, esta aparente paradoja tiene su origen en la difracción 





que ocurre en los bordes de esfera [BH98] y es un resultado conocido. 
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Figura 3.1: Gráficos correspondientes a la Figura 2 de |XF'S14] donde se representa la sección 
eficaz de dispersión dividida por la sección eficaz geométrica para una esfera de material 
dieléctrico con distintos valores del índice de refracción en función del factor de tamaño kk. 
Cálculo realizado con el programa desarrollado para este trabajo 


e Esfera dieléctrica con recubrimiento de grafeno [IC M15]: 
Con el fin de comparar con algunos de los resultados del artículo [P'CM15] para una es- 


fera de material dieléctrico recubierta de grafeno dentro del régimen cuasi-estacionario 





(longitud de onda mucho mayor que las dimensiones de la esfera), derivamos una ex- 





presión asintótica del coeficiente de transmisión del modo magnético, TF, utilizando 





las expresiones asintóticas de las funciones de Bessel para argumento pequeño. 


3.1. VERIFICACIÓN DEL FUNCIONAMIENTO DEL CÓDIGO 


Relative scattering section » 


20 23 


[ej 
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o 
al 
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_ 
al 


Relative scattering section 


(2) 


Relative scattering section 
Oo — N (69) 
SS 
l 
N 
oO 
! 


43 


(3.1) 


(3.2) 
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qe ¡(22 — 1)! 


hee) == on A 


(3.3) 


Vamos a utilizar Ra para el radio de esta esfera, de manera que e; = €2 y gi¡(w) =0. 


Entonces, el coeficiente A72, definido en (2.127) resulta cero y se recupera la forma del 





coeficiente de transmisión de Mie pero con una capa conductora, de la cual da cuenta 


ga(w). 


Para una capa de grafeno de radio Ka en el límite cuasi-estacionario se obtiene: 


(+1) [es +€ + ((+ Doslw)] kr RA 


(3.4) 
A partir de esta expresión asintótica se puede calcular la polarizabilidad multipolar 


[TCM15| dada por: 
¿AE OME ss 


Se tiene entonces: 


pa Mesta ia Jo (6) 


EA Dl ON (3.6) 


En el régimen cuasi-estático, la condición de resonancia de los modos plasmónicos 
viene dada por los polos de «1. Se compara el resultado obtenido en la bibliografía 
utilizando la aproximación Uegi Y keapImlar) + (6pi)*k4,.Ja1|? [TCM15) con el re- 
sultado obtenido usando el código desarrollado. Las figuras 3.3 y 3.4 muestran las 
eficiencias de extinción obtenidas en el presente trabajo y en el artículo mencionado 


respectivamente. 


Esfera dieléctrica con recubrimiento de grafeno [BYZ15]|: 

En el artículo [BYZ15] se analiza la absorción de ondas electromagnéticas por esferas 
de diferentes materiales dieléctricos recubiertas de grafeno. Allí se estudia absorción 
no sólo para una esfera sino para dos esferas separadas y para conjuntos de esferas. Pa- 
ra verificar el funcionamiento del código desarrollado para este trabajo, reprodujimos 
los resultados correspondientes a la eficiencia de absorción de una esfera dieléctrica 


recubierta con grafeno, cuando se varían diferentes parametros. 
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Figura 3.3: Sección eficaz de extinción divida por la sección eficaz geométrica para una capa 
esférica de grafeno de diferentes radios, correspondientes a los de la Fig.3 de [TCM15], con 
medio interior de índices n=1 (a) y n=2,4 (b). Calculado utilizando el código desarrollado. 





El código Python desarrollado calcula las secciones eficaces y eficiencias de dispersión 
y extinción, a partir de las expresiones (2.148) y (2.152), calculadas utilizando los coe- 


ficientes de transmisión, y luego se obtiene la potencia absorbida utilizando (2.141). 








Los resultados obtenidos utilizando el programa desarrollado se pueden ver en la 3.5, 


estos resultados fueron comparados con los obtenidos en [BYZ15) que utiliza una apro- 





ximación de temperatura cero para bajas temperaturas que se pueden ver en la figura 


3.6. 


Luego de haber demostrado que los resultados obtenidos con las herramientas desarro- 
lladas en esta tesis concuerdan de manera aceptable con los resultados disponibles en la 


literatura, en el captulo siguiente obtendremos algunos resultados novedosos para arreglos 
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Tex / TRE 





0.4 0.6 
hw [eV] 


0 0.2 


Figura 3.4: Tomado de [I'CM15], se corresponde con el cálculo de la figura 3.3 utilizando la 


aproximación cuasiestacionaria mencionada. 


utilizando una y dos capas de grafenos concéntricas. 
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R= 100nm 
R=1yum 
R=5um 
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Figura 3.5: Reproducción de los resultados mostrados en la figura 2 de [BYZ15]. En la 
subfigura (a) se tiene la eficiencia de absorción para una esfera de constante dieléctrica e = 
3,9 recubierta con grafeno de potencial químico y =1,0eV y una frecuencia de colisiones 
Ye = 1meV para distintos radius de la esfera. En (b) se utiliza un radio de 100 nm, y. 
lmeV, e = 3,9 y se varía el potencial químico. La subfigura (c) muestra los resultados 
variando la constante dieléctrica para R = 100nm, ig = 1eV y y. = 1meV. En (d) se 
puede ver cómo varía la absorción cuando cambia la frecuencia de colisiones, dejando fijo el 
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radio en 100 nm, el potencial químico en 1,0eV y la constante dieléctrica en 3,9 
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Figura 3.6: Figura tomada de [BYZ15] con los resultados obtenidos por Yang et al. corres- 
pondientes a la figura 3.5 


Capítulo 4 


Resultados 


En las secciones 2.5 y 2.6 hemos hallado las expresiones de los campos electromagnéticos 
dispersados por dos esferas concéntricas de grafeno sobre las que incide una onda plana 
monocromática. Estas expresiones son funciones de la frecuencia incidente, de la geometría 
de la configuración, y de las características constitutivas de los medios así como de las 
superficies de grafeno. 

En esta sección presentamos resultados numéricos relacionados con las características 
esenciales de los modos plasmónicos en los sistemas núcleo-recubrimiento (core-shell). Estos 
sistemas resultan de interés debido a que ofrecen nuevos grados de libertad para sintonizar 
las resonancias plasmónicas que surgen de la hibridización de modos plasmónicos en una 
sola capa. Por este motivo, este capitulo comienza con resultados obtenidos con una sola 


capa de grafeno. 


4.1. Resultados utilizando una capa esférica 


En la figura 4.1 se grafica la sección eficaz de dispersión para una esfera metálica de radio 
R = 50 nm recubierta de grafeno e inmersa en vacío sobre la que incide una onda plana, en 
función de la frecuencia incidente. 


La permitividad eléctrica interior e£¡(w) está dada por la ecuación (4.1) con £x% = 1, 





Vm = 0,01eV y diferentes valores de frecuencia de plasma (w, : 0,5eV,0,7eV y 0,9eV), 
mientras que el valor de la conductividad superficial 0(w) correspondiente al grafeno se 


obtuvo de la expresión de Kubo (1.7) con T' = 300K, y. = 0,1 meV y 1. =0,5eV. 
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w? 


e(w) = € = — E (4.1) 
w(w+ yD) 

Se muestra también la sección eficaz de dispersión para las mismas esferas metálicas 
pero desnudas, así como para una esfera de grafeno vacía en su interior (curva negra, donde 
€1 == =H=1) 

Las frecuencias de excitación están en el rango entre 24 'THz (longitud de onda incidente 


12400 nm) y 170 THz (1770nm). Las curvas de la figura muestran que con grafeno las reso- 


nancias se corren hacia el azul en comparación con las resonancias de las esferas desnudas. 





Esto es consistente con el hecho de que el recubrimiento de grafeno induce un aumento en 


la densidad de portadores de carga en la superficie metálica. 


Solo 
grafeno 


Odispersion 
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Figura 4.1: Frecuencias de resonancia plasmónica para una capa esférica de grafeno y esferas 
metálicas de distinto valor de frecuencia del plasma con y sin recubrimiento de grafeno. 


La figura 4.2 presenta la sección eficaz de dispersión para una esfera metálica de radio 
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R = 50nm recubierta de grafeno (w, = 0,5eV) e inmersa en vacío sobre la que incide 
una onda plana, para diferentes valores del potencial químico (1. : 0,3eV,0,5eV,0,8eV y 
1,1eV). Esto ilustra la capacidad de ajuste de las resonancias, demostrando que estas se 
pueden controlar con la ayuda de un voltaje de compuerta, el cual determina el valor del 


potencial químico. 


dispersión 





Figura 4.2: Sección eficaz de dispersión dividida por la sección eficaz geométrica para una 
capa esférica de grafeno en aire y diferentes valores del potencial químico pe. 


La figura 4.3 muestra las secciones eficaces de extinción correspondientes a una esfera 


de 50 nm de radio y a otra de 200 nm; ambas de permitividad e; = 1 y recubiertas de gra- 





feno (con iy = 0,9eV). Como es de esperarse, debido al aumento de las dimensiones de 
las partícula, la frecuencia de resonancia se corre a frecuencias menores para una esfera de 


radio mayor. 


Se puede apreciar que cuanto más parecidos son los radios de las esferas individuales, 


más cercanas son las curvas de resonancia, y esto es lo que se observa claramente en la figura 
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Figura 4.3: Sección eficaz de extinción correspondiente a una esfera recubierta de grafeno 
de radio 50 nm (verde), otra de 200 nm (azul), en ambos casos con permitividades interior 
y exterior igual a uno y conductividad dada por un potencial químico uu. = 0,9eV. 


4.4. En la misma se presentan las curvas de extinción como las de la figura 4.3, pero esta 
vez para esferas de radios parecidos (una de radio 198nm y la otra de 200 nm). Observar 


que las curvas son practicamente coincidentes. 





Finalmente, se realizaron gráficas del módulo del campo eléctrico sobre la superficie para 
una superficie esférica de grafeno autosostenida, esto es, con aire o vacío tanto en el espacio 
interior como exterior a la superficie. La eficiencia de extinción para radio R = 1000 nm, 


con un potencial químico 4. = 1,3eV, una frecuencia de colisiones y. = 0,0001eV y a una 





temperatura de 300 K puede verse en la figura 4.5. Allí aparecen una serie de picos que 





corresponden a resonancias de diferentes modos multipolares, donde el orden del multipolo 


crece a medida que aumenta la frecuencia de excitación y se aprecia que la intensidad de 
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Figura 4.4: Sección eficaz de extinción correspondiente a una esfera recubierta de grafeno 
de radio 198nm (verde), otra de 200 nm (azul) con conductividad dada por un potencial 
químico 1. =0,9eV, las permitividades en los espacios interior y exterior son iguales a 1. 


los picos decrece a medida que aumenta el orden. 


La figura 4.6 muestra el módulo del campo eléctrico a un radio 1,01R donde R es el 


radio de la esfera, esto es, en el espacio exterior a la superficie, justo por encima de ella. 





La incidencia es en la dirección Z y el campo eléctrico está polarizado en la dirección X. 
Se distinguen en los diferentes gráficos los polos (en rojo) donde el campo es más intenso. 
Debido a que la intesidad de la resonancia disminuye con el orden multipolar, a medida que 


aumenta este orden se vuelve más difícil identificar los polos, en el caso octupolar algunos 





polos no llegan a resolverse. 
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Figura 4.5: Eficiencia de extinción para una capa esférica de grafeno autosostenida. Se 
pueden apreciar los picos correspondientes a las excitaciones de distintos multipolos. El radio 
de la esfera es 1 ym y para su conductividad se utilizó un potencial químico 4. = 1,3eV, 
una frecuencia de colisiones y. = 0,0001eV y una temperatura de 300 K 


4.2. Resultados utilizando dos capas esféricas concéntri- 


cas 


Para ver la interacción entre dos capas de grafeno se presentan en la figura 4.7 las curvas 
de la figura 4.3 junto con la curva correspondiente a la sección eficaz de extinción para dos 
superficies esféricas juntas. En este caso se observa claramente que los picos de resonancia 
dados por las dos cáscaras de grafeno concéntricas juntas corresponden a la unión de los 
picos de ambas esferas individuales. 

Sin embargo, cuando se observa la interacción entre dos capas de grafeno muy cercanas 
(como es el caso de los radios de la figura 4.4) el efecto difiere del presentado en la figura 4.9 
para radios más alejados. Aquí la curva para la sección eficaz de extinción de dos cáscaras 
de grafeno de radios parecidos está claramente corrida respecto a las curvas individuales 
(ver figura 4.8). 

La interpretación de lo observado en la figura 4.8 es que al estar muy juntas ambas 
superficies de grafeno, el efecto en la sección eficaz es similar al que produce una única su- 


perficie, pero con una conductividad mayor a la de las dos capas individuales. Este resultado 


4.2. RESULTADOS UTILIZANDO DOS CAPAS ESFÉRICAS CONCENTRICAS 313) 
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Figura 4.6: Intensidad del campo eléctrico a un radio de 1,01 veces el radio de la esfera. 
El gráfico de cada resonancia aparece en dos orientaciones diferentes. En el sistema de 
ejes, el eje azul representa la dirección Z a lo largo de la cual incide la onda plana, el eje 
rojo representa la dirección X que contiene la polarización del campo y el eje verde es la 
dirección Y. Se puede ver que en la resonancia dipolar el campo tiene dos polos que son 
máximos de intensidad, en el caso cuadrupolar cuatro, etc. En la resonancia octupolar los 
polos se distinguen con mayor dificultad debido a que la intensidad de la excitación es mucho 
menor que en los otros casos y que la presencia de la cola de otras resonancia desdibuja la 
distribución. 
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Figura 4.7: Sección eficaz de extinción correspondiente a una esfera recubierta de grafeno 
de radio 50nm (verde), otra de 200nm (azul) con conductividad dada por un potencial 
químico 4. = 0,9eV, y sección eficaz de extinción para ambas esferas juntas (rojo) (con 
permitividades en los medios £, = €2 =€3 = 1). 


es consistente con el análisis de dos capas que se hace en la sección III de ([DSO13]). 
Para constatar esto se presenta en la figura 4.9 la curva de la sección eficaz de las mismas 
dos esferas juntas de la figura 4.8, las que tienen radios 198 nm y 200 nm, ambas de grafeno 


con conductividad dada por un iy = 0,9eV, junto con la curva correspondiente a una 





única esfera (de radio 200 nm), pero en este caso con conductividad correspondiente a un 
Mg = 1,8eV. 
En la figura 4.10 se presentan curvas de extinción para dos esferas concéntricas, de 


grafeno, de Ri = 190nm y Ri = 200nm, ambas con ig = 0,9eV, con aire en la región 





interior y exterior (e; = £3 = 1), pero para varios valores de permeabilidad en la zona 


intermedia, donde se dejó fija la parte real de la misma (Re(e2) = 1,5) y se varió la parte 
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Figura 4.8: Sección eficaz de extinción correspondiente a una esfera recubierta de grafeno de 
radio 198 nm (verde), otra de 200 nm (azul) con conductividad dada por un potencial químico 
pic = 0,9 eV, y sección eficaz de extinción para ambas esferas juntas (con permitividades en 
los medios e1 =€2 =€3 = 1) (rojo). 


imaginaria (Im(e2) : 0, 0,01, 0,05 y 0,15). 
La parte imaginaria de e2 da la pérdida de energía por absorción en la correspondiente 
zona; esto se observa en que al aumentar Im(e2) se van ensanchando los picos de resonancia, 


pero sin cambiar la posición de los mismos. 


4.3. Estudio de los campos cercanos 


Para complementar las curvas de extinción que se dan para dos cáscaras esféricas inter- 
actuando, se presentan en la figura 4.12 los campos cercanos (módulo de campo eléctrico) 


correspondientes a dos esferas de grafeno de radios Ri = 0,9 ym y Ra = 1 um, ambas con 
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Figura 4.9: Sección eficaz de extinción correspondiente a dos esferas recubiertas de grafeno 
juntas de radios 198 nm and 200 nm (rojo) con conductividad dada por un potencial químico 
pi. = 0,9eV (con permitividades en los medios e, = €2 = €3 = 1) y sección eficaz de extinción 
para una única esfera de radio 200 nm con potencial químico u. =1,8eV (azul). 


Ma = 1eV, e inmersas en aire (e, = €2 = €3 = 1). En las figuras presentadas la onda 
excitadora incide desde abajo y en cada caso la excitación es a la frecuencia de resonancia 
correspondiente al pico indicado en la figura 4.11 (5,85 THz (a), 9,89 THz (b), 19,9 T'Hz 
(c), 26,95 THz (d), 31,76 THz (e) y 35,5 THz (f)). En los gráficos se incorporaron versores 


indicando la componente del campo paralela al plano del corte. 


Claramente las simetrías que presentan las figuras (a) y (c) corresponden a excitaciones 





dipolares, (b) y (d) a cuadrupolares, mientras que (e) es hexapolar. Hagamos notar que las 
dipolares presentan un campo uniforme en el interior. Por otra parte, (a) y (b) presentan 


sus valores máximos de intensidad en la zona intermedia mientras que (c) y (d) lo hacen en 





la zona externa. El caso de (e) es mas complejo aunque puede decirse que la mayor parte 
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Figura 4.10: Eficiencia de extinción para una partícula con núcleo de aire o vacío y dos capas 
de grafeno separadas por un medio con permitividad compleja. En las diferentes curvas se 
varía la parte imaginaria de la permitividad. Se puede apreciar que el aumento de esta parte 
imaginaria, que al ser positiva representa absorción en el medio, tiene el efecto de ensanchar 
los los picos de resonancia sin alterar sus posiciones. 





del máximo de intensidad se da en la zona interna a la esfera mas chica. 
Por último, observar que la escala de intensidad no es la misma en cada figura, disminuyendo 
el valor máximo al aumentar el multipolo, y concordando aproximadamente con la altura 


relativa entre picos de la figura 4.11. 


De forma similar a la anterior, en la figura 4.13 se muestran los campos cercanos corres- 
pondientes a los picos dipolares presentados en la figura 4.9. El caso (a) es para una única 
esfera de radio R = 200 nm, de grafeno, donde la conductividad viene dada por un potencial 
químico y = 1,8eV, y el caso (b) es para dos esferas de radios R¡ = 198 nm y Ra = 200 nm, 


donde el grafeno tiene en ambas capas una conductividad dada por y = 0,9eV. 
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Figura 4.11: Eficiencia de extinción para dos superficies esféricas de grafeno de radios R¡ = 
0,9 um y Ra = 1 ¡um con un potencial químico 4, = 1eV y e, = €2 = €z = 1. Los picos 
que pueden verse, (a), (b), (c), (d), (e) y (f), representan resonancias correspondientes a 
diferentes excitaciones multipolares. 





Es clara la similitud en ambos casos, aunque hay que hacer notar que la escala de valores 
de intensidad no es la misma. 
La figura 4.14 muestra la sección eficaz de dispersión en función de la frecuencia de la 


onda incidente para dos esferas recubiertas de grafeno de radios R¡ = 195 nm y Raz = 200 nm, 





para cuatro valores diferentes de potencial químico en las capas de grafeno. En todos los 





casos los tres medios adyacentes son aire. Se puede observar que la disminución del potencial 


químico implica un corrimiento al rojo de la frecuencia de resonancia. 
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Figura 4.12: Mapa de campo eléctrico cercano (|E|) para los cinco multipolos indicados 
en la figura 4.11 excitados desde abajo con una onda plana monocromática, de frecuencia 
correspondiente al modo excitado. Las superficies esféricas se indican en verde. Se incluyen 
vectores indicando la componente del campo contenida en el plano (y,z). En la versión 
electrónica se sugiere ampliar la figura para verla con mayor detalle. 
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Figura 4.13: Intensidad de campo eléctrico cercano para: (a) Una capa esférica de radio 
de 200 nm utilizando un potencial químico de 1,8eV y (b) Dos capas esféricas de radios 
de 198nm y 200 nm utilizando un potencial químico de0,9eV. En la versión electrónica se 
sugiere ampliar la figura para verla con mayor detalle. 
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Figura 4.14: Sección eficaz de extinción correspondiente a dos esferas recubiertas de grafeno 
juntas de radios 195 nm y 200 nm (con permitividades en los medios £1 = €2 = €3 = 1) para 
diferentes valores de potencial químico (el potencial químico en ambas capas siempre es el 
mismo). 
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Capítulo 5 


Conclusiones 


En este trabajo, enmarcado en el estudio de posibles aplicaciones de las estructuras 
core-shell con grafeno en el rango de los terahertz, estudiamos la excitación de plasmones en 
estas estructuras. Mientras que en el caso metálico, en el cual las paredes son impenetrables 
en el mencionado rango de frecuencias, y en el caso dieléctrico, en el cual las propiedades 
dinámicas las determinan la diferencia de índice entre los medios; en el presente caso las 
propiedades de los plasmones excitados dependen también de las propiedades constitutivas 


del grafeno, dadas por su conductividad. 





Aplicamos las ecuaciones del electromagnetismo macroscópico y a partir de las condi- 





ciones de contorno hallamos de forma rigurosa los coeficientes de transmisión que conectan 
los campos en las distintas regiones dieléctricas separadas por esferas concéntricas de gra- 


feno; en estos coeficientes aparecen términos que contienen explícitamente la conductividad 





del mismo. Con estos coeficientes se hallaron expresiones exactas para los campos electro- 


magnéticos en todo punto del espacio en función de la amplitud de la onda incidente. 


Hemos obtenido resultados numéricos y como control de los códigos numéricos desarro- 





llados hemos verificado que los resultados reproducen los casos particulares presentes en 
bibliografía. A partir de allí, hemos ilustrado aspectos dinámicos de los plasmones que se 


excitan en las dos superficies. 


Aunque el formalismo desarrollado en este trabajo permite tratar situaciones con paráme- 
tros constitutivos generales para los medios interior, intermedio y exterior a las superficies 


esféricas de grafeno, hemos considerado unas pocas situaciones ilustrativas concentrándonos 
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en la interacción entre ambas esferas al variar el radio de las mismas. 
Ilustramos con ejemplos la sintonizabilidad de la frecuencia de resonancia a través del 


potencial químico del grafeno y también el comportamiento de la frecuencia de resonancia 





plasmónica al variar el radio de las esferas y vimos que se presenta un fenómeno de hibridi- 
zación al acercarse los radios de las mismas. Hemos mostrado el confinamiento transversal 


del campo electromagnético asociado a una excitación plasmónica sobre las paredes de una 





esfera. Se observó la absorción dada en los medios, dependiente de la parte imaginaria de la 
permitividad eléctrica. También hemos determinado las características polares de los modos 
plasmónicos £ = 1,2 y 3 mostrando distribuciones de campo cercano. 

Como continuación natural del presente trabajo pueden mencionarse el estudio de la 
dispersión de ondas electromagnéticas por capas esféricas de grafeno cuya conductividad 
no sea homogénea o el estudio de la dispersión por capas planas constituidas de pequeñas 
esferas dieléctricas recubiertas de grafeno. "También resulta de interés estudiar la dispersión 
por esferas constituidas de un número apreciable de capas dieléctricas o metálicas separadas 


por superficies de grafeno. 
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